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Abstract 



<^ ' A nonpolynomial one-dimensional quantum potential in the form of an isotonic oscillator (har- 

cn , monic oscillator with a centripetal barrier) is studied. We provide the non-relativistic bound state 

energy spectrum En and the wave functions ipni^) ™ terms of the associated Laguerre polynomi- 

O i' als in the framework of the Nikiforov-Uvarov method. Under the spin and pseudospin symmetric 

"^ ■ limits, the analytic eigenvalues and the corresponding two-component upper- and lower-spinors of 

the Dirac particle are obtained, in closed form. 

Keywords: Schrodinger equation, Dirac equation, spin and pseudospin symmetry, harmonic 
K^ ■ oscillator, isotonic oscillator, self-adjoint operator, Nikiforov-Uvarov method. 

m 



PACS numbers: 03.65.Pm; 03.65. Ge; 02.30. Gp 



* E-mail: lsikhdair@neu.eduTtrl 
^E-mail: lsever@metu.edu. trl 



I. INTRODUCTION 

The exact solutions of the nonrelativistic and relativistic wave equations are only possible 
in a few simple potential cases such as the Coulomb, the harmonic oscillator, the pseudo- 
harmonic, isotonic potentials and others [1-4]. The most interesting and best known system 
inside this small family is the harmonic oscillator whose energy spectrum consists of an 
infinite set of equidistant energy levels. Many other oscillators, as for example harmonic 
oscillators perturbed by a term containing a fourth or a sixth power in the coordinate, have 
been extensively studied. Nevertheless, it is known the existence of another solvable one- 
dimensional model which shares several interesting properties with the harmonic oscillator, 
the so-called isotonic oscillator whose spectrum coincides with that of the harmonic oscil- 
lator [5] . This potential model is considered as one of the most used models for the study 
of the dynamics of nonlinear systems. In particular, the energy spectrum of the singular 
potential is considerd as an isomorphous to the harmonic oscillator spectrum. Very recently. 
Fellows and Smith [6] have used the ideas of the factorization and supersymmetric (SUSY) 
approach [7,8] to study the singular superpotentials and they proved that most of these 
oscillatory potentials are other partner potentials of the harmonic oscillators and derived 
an infinite set of exactly soluble potentials. A discussion of the supersymmetric connection 
between harmonic and isotonic oscillators can be traced in Ref. [9]. 

The isotonic potential takes the form. 



ImwV + H 

2 2a;2 



Uisotix) = Uo{x) + Ug{x) = -MwV + -^, X ^ 0, g>0, (1) 



where u is the angular frequency of oscillator, Mu"^ = K in classical mechanics and g = 
m{m -\- 1). The common feature of this potential is that it consists of a harmonic term plus 
an additional rational function (centripetal barrier) which falls off at infinity like a constant 
g times 1/x'^ with one regular singularity at x = along the whole domain — oo < a; < oo.[l] 
It is clear from (1) that Ug{x) exhibits strong singularity when a; = 0, so that the wave 
functions must vanish at such a point. It is worthy to note that the Hilbert space associated 
to the models with Ug{x) is narrowed as compared to the Hilbert space of the hamiltonian 
with potential Uq{x). In addition, the domain of the hamiltonian associated to harmonic 



[1] Centrifugal barrier does not make physical sense in one-dimension. It is often used for such singular terms 
in a potential. 



oscillator extends along the whole real axis (—00,00), however, for potential (1) reduces 
itself to the half-line (x > or x < 0). The harmonic oscillator Uo{x) and the isotonic 
oscillator Ujsot{x) are plotted in Figure 1 for the sake of comparison. The aim is to show 
that both curves coincide for wide range of a; > except for values of x in the neighborhood 
of the vertical asymptotic line x = where the isotonic oscillator goes to infinity. They 
have also identical spectrum, however, shifted by two units. We consider this potential in 
the interval (0, 00) as in the case of (1). Authors of Refs. [7-9] have solved the Schrodinger 
equation for the potential (1) with u = M = 1 and g = 2 {m = 1) using SUSY approach 
[7,8] on the assumption that the superpotential W{x) obtained from the wave function of 
the harmonic oscillator 0i(x) can generate twice the potential (1) (see e.g., [6]). The singular 
term in the above potential is often called the centripetal barrier potential. However, the 
centripetal barrier potential here makes no physical sense in one dimension since the term 
m{m + l)/x'^ singularities are often related to the radial equation for the three-dimensional 
harmonic oscillator. 

In this article we set up to present a study of the exact analytic nonrelativistic bound 
state energy spectrum and the correspoding wave functions in terms of the associated La- 
guerre polynomials L'^{z) (or the Kummer confiuent hypergeometric function, M{a,b,z) = 
iFi {a;b;z)) by applying the Nikiforov-Uvarov (NU) method [10]. Overmore, we extend 
our study to investigate this potential model in the context of the spin and pseudospin 
symmetric Dirac equation [11-17]. In the presence of the spin symmetry S ^ V = Uisot{x) 
and pseudospin symmetry S ~ —V = Uisot{x), we investigate the exact s-wave bound state 
energy eigenvalues and corresponding upper and lower spinor wave functions in a systematic 
form [12-15]. We also show that the spin (pseudospin) symmetric Dirac solutions can be re- 
duced to the S = V = Uisot{x) {S = —V = Uisot{x)) in the presence of exact spin symmetry 
A = (pseudospin symmetry S = 0) limitation [16,17]. Overmore, the solution of the Dirac 
equation can be easily reduced to it's nonrelativistic limit if one applies an appropriate map 
of parameters. 

The rest of the paper is organized as follows. In section 2, we apply the NU method 
to solving the Schrodinger and Dirac equations with an exactly solvable isotonic oscillator 
to obtain the eigenvalues and eigenf unctions in a systematical way. We also compare the 
non-relativistic solution with the existing one obtained by applying the SUSY approach. In 
section 3, we make our summary and conclusions. 



II. BOUND STATE SOLUTIONS 



A. Schrodinger case 



We start with the one-dimensional single-particle Schrodinger equation [18,19]: 



Hisoti'ni'^x) 



+ f//sot(±a;) 



^^{±X) = En^^{±x), 



(2) 



2Md{±xy 

The potential in Eq. (1) is invariant with respect to the inversion, i.e., Ujsot{~x) — )■ Ujsot{x) 
as well as the hamiltonian Hjsot- Therefore, the Schrodinger equation should have even and 
odd solutions. The domain of the harmonic oscillator hamiltonian H extends along the 
whole real axis — oo < x < oo, however, the partner Hjsot exhibits a strong singularity 
at the origin. The space breaks up into two disjoint regions (a; > or a; < 0) without 
communication between them since the wave functions vanish at the regular singularity 
X = (i.e., ipni^) = 0) and at the irregular singularities ±oo (i.e., ■?/'„(±oo) = 0). In this 
respect, we should restrict the hamiltonian to the interval (0,oo), this is exactly the same 
situation that occurs in isotonic potential. In addition, if we set x — )■ ix, then (2) becomes 



2Mdx^ 



+ Uisot{x) 



^„(ix) 



-En^ni'^x) 



which is the original equation with the irrelevant change in the eigenvalues En — )■ —En- If we 
set a; — )■ ix in the isotonic wave function, we get perfectly good ipni'i-x) which can generate 
the superpotential W{x) and other hamiltonians [6]. 

To solve Eq. (2) by NU method, we perform a straightforward algebra to reduce it into 
the following simple form: 



<(^) + 






with 



2MEn 



(3 



Mu 



a 



Mg 



h^ 



> 0. 



Let us restrict ourselves to the positive half-line {x > 0) and in terms of new variable s 
(0 < s < oo), we obtain 



(3) 



(4) 



X 



(5) 



where we have used ipni^) = '?/'„(s). Now, if we compare the above equation with the following 
generalized hypergeometric-type equation with a parametrization of real variables s = s{x): 

<{s) + ^-^Ms) + ^V^„(s) = 0, (6) 

where 

^n(s) = l^(s)y„(s), (7) 

and where a{s) and a(s) are two polynomials, at most of second-degree, and t{s) is at most 
of first-degree polynomial, then it follows that: 

t{s) = 1, a{s) = 2s, a{s) = -(5'^s^ + £„s - a. (8) 

To apply the NU method [10,20], we calculate the function 7r(s) defined by 



a'is] -T[s] 



TT S = ^ ^ ± 



a [s] — Tis] 



2 



a{s) + ka{s) 



i ( 1 ± \JaP^s^ + 4 (2A; - £„) s + 4a + 1 j 



(9) 
and also seek for a physical value of k that makes the discriminant of the expression under 



square root, in the last equation, to become zero (i.e., 2A; = £„ ± /3y/T+4a, a > —1/4). 
Hence, there is no bound solution in the region a < —1/4. The model becomes unphysi- 
cal in the region (— oo,— 1/4), since the spectrum is not bounded from below. Upon the 
substitution of the value of k into the above equation, we obtain the following suitable 
solutions: 

7r(s) = 2 (l + v^lTl^) - /3s, (10) 

and 

k = ^(^en-PVmaj. (11) 

With regard to Eqs. (8) and (10), we can calculate the function t{s) = t{s) + 27r(s), taking 
into consideration the bound state condition which has to be established when t'(z) < 0, as 



t{s) = 2 + VT+Ia - 2/3s and t'{s) = -2(3 < 0. (12) 

According to the method, in order to find the energy equation from which one calculates 
the energy eigenvalues, we need to find the values of the parameters: A = k + tt'{s) and 
A = A„ = —nT'{s) — ^n{n — 1) (y"{s), n = 0, 1, 2, ■ ■ ■ , as 

X = -(en- l3Vl + 4a) - l3 and A„ = 2n/3, n = 0, 1, 2, ■ ■ ■ . (13) 



Using the relation A = A„ and the definitions of parameters in Eq. (4), one finds that the 
energy eigenvalues of the isotonic oscillator are 



Er,,g = hw\2n+1 + ^^1 + ^Y n = 0,l,2,--- , (14) 

which is identical to the results of Ref. [21] (see p. 3). Thus, we have exactly solved 
the isotonic Hamiltonian (2). In the limit that (? — )■ 0, the relation (14) reduces to En = 
hu! (2n + I) , n = 0, 1, 2, ■ ■ ■ which is identical to the s-wave solution of the three-dimensional 
Schrodinger equation with harmonic oscillator potential (cf. Eq. (34) when / = 0). 

Overmore, using the conventions of Ref. [9] (cf. Eqs. (41) and (43) therein), we may 
take g = 2 (i.e., m = 1) for easy of notation, the spectrum of Eq. (1) {in h = M = u = 1 
units) reads as 

E„ = 2n+| n = 0,l,2,---. (15) 

It is noticed that energy spectrum in the previous equation is half the spectrum of Eq. (41) 
in Ref. [9] (see Eq. (43) in [9]). The odd solutions under the inversion x — )■ —x (the negative 
half-line, x < 0) has same energy spectrum as the even ones given in Eq. (14) due to the 
invariance of isotonic potential under this inversion. Hence, the energy relation in Eq. (14) 
holds for the whole domain — oo < a; < oo. 

Thus, all the other eigenenergies are given by 

E2n = Eo + 2nu, n = 0, 1, 2, ■ ■ ■ , (16) 

and the energy spectrum is equidistant since 

E2n+2 = E2n + 2u. (17) 

Nevertheless, the height Ai? = 2oo of the energy steps is twice that of the simple harmonic 
oscillator Uq. In fact, it seems as if half of of the states (those with an odd number of nodes) 
have disappeared. 

Let us now turn to the calculations of the normalized wave function. The first part of 
the wave function in Eq. (7) is found through the relation [10,20]: 

n{s) = exp( f 44rfs) = s^+^^ exp (-^/3s) , (18) 



with 



^ = lv/rT4^>l, (19) 



and the calculation of the weight function is found through the relation, 

Pis) = ^ exp (J 44^5 J = s« exp (-/3s) , (20) 

leading to the calculation of the other part of the wave function; namely, |/„(s) which is 
a hypergeometric type function whose polynomial solutions are given by the Rodrigues 
relation: 

Vnis) = Anp~\s)-^ Kis)pis)] = s'^ sxp (/3s) ^ (s"+« cxp (-/3s)) = L«)(/3s), (21) 

:where L'^{y) is the associated Laguerre polynomials. Therefore, the even solution of the 
wave function satisfying Eq. (7) is [22,23] 



It should be noticed that the change oi x ^ ix in the eigenvalue equation (2) for the isotonic 
oscillator changes eigenvalues En — > —En. Then, if ipni^) is the eigenfunction corresponding 
to the eigenvalue En, then the eigenfunction ■?/'„(ix) will be normalizable only if —En is 
in the point spectrum of this Hamiltonian. However, the ■?/'„(ix) would not be good wave 
functions if their spectrum —En lies not within the Hamiltonian range, then ■?/'„(ix) is not 
normalizable eigenfunction. At first, it behaves like exp Q/Ja;^) at large x but this is not 
relevant here as one writes the explicit form of Ln {(3x'^) for even n. On the other hand, the 
odd solutions for which the part of wave function corresponding to —x have opposite signs 
and exist as 

iJni-x) = Mn (-a:)^-«exp (-\px'\ Lf{[5x^). (23) 

The two linearly independent solutions (wave functions) given by Eqs. (22) and (23) for the 
even and odd solutions, respectively, need to be normalizable in the range (0, oo). However, 
the odd solution is not normalizable in the region (0, oo) as one can see in Eq. (23). Indeed, 
the operator in (3) is not essentially self-adjoint for —1/4 < a < 3/4 and its most general 
square-integrable solution behaves near the singularity as a linear combination of x^^"^^^ and 
^1/2-5 (See Eqs. (22) and (23)). [1] In this respect, we clarify this point by analyzing the 
behaviour of the isotonic potential (1) in terms of the parameter a [24,25], three different 
regions appear, namely. 



[1] We would like to thank one of the referees for drawing our attention to this point. 



• In the range a G (—00, —1/4) the model becomes unphysical since the spectrum is not 
bounded from below (see Figure 1) [24]. 

• When a G ( — 1/4, 3/4), the singularity is not strong enough to make the wave functions 
(22) and (23) vanish at a; = 0. Indeed in this region both linearly independent solutions 
(wave functions) are normalizable since ^ = and 1. This is the reason why it 
is necessary to select, from the continuous family of self-adjoint extensions by the 
differential operator, the self-adjoint extensions which correctly describes the physical 
system under consideration [24]. The wave functions pass across the singularity point 
a; = and the model extends itself again along the entire region, that is; (—00, 00) . 

• Physically, if we consider the range a G (3/4, 00), the singularity acts as an impentrable 
barrier, thus dividing the space into two independent regions, that is; a; < and 
a; > 0. The wave functions must vanish at a; = which provides an absolute lack 
of communication between the two regions of space (i.e., the negative and positive 
half-lines) and the wave functions (22) and (23) in this case become normalizable 
[24,25]. 

Alternatively, notice the Laguerre polynomial can be expressed in terms of the Kummer 
confluent hypergeometric functions as [23] 

^"n{z)=^^^,F, (-n;p + l;z)), (24) 

where 

/ , X -, « a(l + a) 2 a(l + a)(2 + a) o ^. ,. ,^. 

,F, (a; 6; .) = 1 + -. + ^1^.^ + 66(1 + t)(2 + 6? + ''^ ' ^^^ 

Using the notations of other authors {h = M = 1)[21] and putting (5 = u and ^ = m -{- 1/2, 
the even wave function solution in Eq. (22) becomes [23] 

^„(x) = Ar„a;^+" exp i --ujx'^ J iFi ( -n; m + -; wa;M , n = 0, 2, 4, ■ ■ ■ , 



r (m + I) V n\ 

which is identical to Eq. (42) in Ref. [9] when we set m = 1 so that the isotonic potential 
given by Eq. (41) in [9] is twice the potential (1) in the present work. On the other hand. 



the odd solution for which the part of wave function corresponding to —x can be expressed 

as 

^^{-x) = (-l)i+™^„(a:), n = 1, 3, 5, ■ ■ ■ . (27) 

Hence, if we take m = 0,2,4, ■■■ (even real integer), we find the wave function ■?/'„(— a;) 
being an odd (antisymmetric) function of x [i.e., ■?/'„(— x) = —iIJj^{x)]. Overmore, if we take 
m = 1, 3, 5, ■ ■ ■ (odd real integer), we find the wave function corresponding to negative values 
of X is identical to the wave function corresponding to positive values, i.e., ■?/'„(— x) being 
an even function (symmetric) with ■?/'„(— a;) = iIJj^{x). In case if m is a real number but not 
integer yielding a complex wave function in the negative half-line which is not normalizable. 
On the other hand, the energy levels of the one-dimensional Schrodinger equation for the 
harmonic oscillator Uo{x) [1]: 

En= (n+^jhu, n = 0,l,2,--- , (28) 



and the well-known wave functions [26] 

1/2 

"I / W 



2"n! V TT 



H^{f3x')exp{~Px'), P = ^, (29) 



where Hn{y) = {—!)"' exp[y'^)j^ exp{—y^) represent Hermite polynomials. 

For further discussions on the isotonic potential, we present the energy eigenvalues and 
the corresponding wave functions of the lowest three states, 

Eo = oo(^ + mY (30a) 

^^(x) = iVoa:^+™ exp ( -^rOoxA , (30b) 



Ei=\J + mj u, (31a) 

^,(.) ^ iV,.'-»exp (-i..^) (1 - pJ^.^ + ^__||__^.^) , (32b) 



and 



respectively and the normalization factors are calculated as 



- £. / i-vi_L.^ //~i . r-i \ / /yy^J-A 



The isotonic ground state wave function, ip^ix) in (30b) is compared with the corresponding 
harmonic oscillator wave function, (/)q{x) in Eq. (29) in Figure 2. Further, the isotonic first 
two excited wave functions, ipi{x) and ip2{x) in Eqs. (31b) and (32b) are compared with the 
corresponding harmonic oscillator wave functions, 0i(a;) and 02 (^) ^^ Eq. (29) in Figures 3 
and 4, respectively. 

On the other hand, the solution of the three-dimensional Schrodinger equation with any 
arbitrary quantum number / (i.e., harmonic oscillator combined with centrifugal barrier 
potential) provides us 

En,i = hw (2n + I + ^ j, n,l = 0,1,2,- ■■ (34) 

and the corresponding wave functions are given by 

i^nAr, e, V.) = y(^y^^^^^-' exp(-l/3r^)L^V2) (^,2) ^^^(^^ ^^^ (35) 

where Ln (/3^^) is the associated Laguerre polynomial, and Yi^mi^, '■p) is the angular part 
of the wave functions. The order n of the polynomial is a non-negative integer. Thus, the 
exact solution of the isotonic oscillator in Eq. (14) in one-dimension is equivalent to the 
solution of the harmonic oscillator f/o(^) combined with the centrifugal barrier potential 
/(/ + l)/r^, r G (0,00), in three-dimensions given in Eq. (34) when we take h = M = 1 
and g = m{m + 1), where (7 is a real number. That is, En,m = ^ (^1 + |) is equivalent to 
En^i = hcu (n2 + I) , where we have defined rii = 2n + m and n2 = 2n + 1 [m -^ I) for which 
the solutions are defined for positive half-line [x G (0, 00) ^ r G (0, 00)]. 



B. Dirac Case 

We start by writting the two radial coupled Dirac equations for the upper and lower (i.e., 
Fn,K{j') and Gn,K{j'), respectively) spinor components [26,27]: 



(^ + 7) FnAr) = {Mc^ + En. - A) Gn,.{r), 



(36a) 
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(-^ - ^) GnAr) = {Mc' - K« + S) F„,«(r), 



(36b) 



where A = V — S and J] = V + S are the difference and sum potentials, respectively and 
are expressed in terms of vector {V) and scalar (S) potentials. In addition, c ~ 137 is the 
velocity of light. 

In the presence of spin symmetry ( i.e., A = Cg), one gets a second-order differential 
equation satisfying the upper-spinor component [17,28-31] 

FLir) - (^^^^ + AI + 7S) K^r) = 0, (37) 

where 

Al = ^ [Mc^ - En.) , 1 = ^ {Mc^ + Kk - Cs) > 0, (38) 

and K (k + 1) = / (/ + 1) , K, = I for k, < and /t = — (/ + 1) for k, > 0. The spin symmetry 
energy eigenvalues depend on n and k, i.e., En. = E{n, k{k, + 1)). For / 7^ 0, the states with 
j = I ± 1/2 are degenerate. Further, the lower-spinor component can be obtained from Eq. 
(36a) as 

where En. 7^ —Mc^, i.e., only real positive energy states exist when Cs = (exact spin 
symmetric case). 

On the other hand, under the pseudospin symmetry ( i.e., S = Cpg), one obtains a 
second-order differential equation satisfying the lower-spinor component, 

G'Lir) - (Hin^ + Al - ^a] G.Ut) = 0, (40) 



^2 ' PS 



where 



Al = 7 {Mc^ + E„,) , 7 = -^ {Mc^ - Enn + Cps) , (41) 



and the upper-spinor component Fn^ir) is obtained from Eq. (36b) as 

where En. 7^ Mc^, i.e., only real negative energy states exist when Cpg = (exact pseudospin 
symmetric case). From the above equations, the energy eigenvalues depend on the quantum 
numbers n and k, and also the pseudo-orbital angular quantum number / according to 
k,{k, — 1) = /(/ + 1), which implies that j = I ±1/2 are degenerate for / 7^ 0. The quantum 

11 



condition for bound states demands the finiteness of the solution at infinity and at the origin 
points, i.e., Fnn{0) = GnniO) = and F„k(oo) = Gn^loo) = 0. 

Let us now study the isotonic potential (1) in the context of spin and pseudospin sym- 
metric Dirac equations. It is well-known that Eqs. (37) and (40) can be solved exactly 
for any k, with the spin-orbit (pseudospin-orbit) centrifugal (pseudo centrifugal) potential 
term. However, we shall study these equations for the s-wave case {k = ±1) for the sake of 
comparison with the nonrelativistic case since m{m + l)/x'^ in the isotonic potential has the 
same behaviour as k,{k, ± l)/r^ in Eqs. (37) and (40). 

1. Spin symmetry limit 

This symmetry arises from the near equality in magnitude of an attractive scalar, S, and 
repulsive vector, V, relativistic mean field, S* ~ V^ in which the nucleon move [12]. Therefore, 
we simply take the sum potential equal to the isotonic potential model, i.e.. 



2 2a;2 



Ui^otix) = -Mco^x^ + -^. (43) 



In the last equation, the choice of S = 2V — ?> Uisot{x) as stated in Ref. [26] allows one to re- 
duce it into its non-relativistic limit under appropriate choice of parameter transformations. 
Further, we take k = — 1 (/ = 0) and in terms of new variable s = x"^ (positive half-plane 
X > 0), Eq. (37) becomes 

K-l{s) + -^fn,-l{s) + ^ [-^'^' - ^l' - P] FnMs) = 0, (44) 

where 



13 = -c/7 and z/ = J-Mu^-f. (45) 

The quantum condition is obtained from the finiteness of the solution at infinity and at 
the origin point. We apply the NU method following the same steps of solution in previous 
section to obtain the expressions: 

t{s) = 1, a{s) = 2s, J{s) = -uh^ - A^s - (3. (46) 

It follows that the functions required by the method for 7r(s), k and t{s) take the suitable 
forms: 

<s) = -z/s + i (l + ^iTa^) , (47) 
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l(^Al + uy^mp], (48) 



and 



t{s) = 2 + VI + 4/3 - 2z/s and t'{s) = -2v < 0, (49) 



respectively, with prime denotes the derivative with respect to s. Also, the parameters A 
and \n take the forms: 

\ = -- (aI + u^l + A[5\ - z/, and A„ = 2nu. (50) 

Using the condition A = A„ followed by simple algebra, we obtain the following transcen- 
dental energy equation. 



(E„,_i - Mc^) ^Mc^ + E„,_i - Cs = haoy/2M ( 2n + 1 + ]:J^ {Mc^ + K,-i - C,) + 1 



\ 2 V fPc^ 

(51) 

where n = 0, 1, 2, 3, ■ ■ ■ and En,-i > Cg — Mc^. One can compute the energy eigenvalues by 

choosing suitable parameters in the symmetric potential. Hence, Eq. (51) shows the energy 

eigenvalues En dependence on n and Cg as well as on the parameters u and M. 

Therefore, using Eq. (51), we compute some energy levels for several values of n (in 
units ^ = c = 1). In the presence spin symmetric limit. Table 1 gives some numerical 
results by taking the following parameters values: M = u = 1.0 fm^^, Cs = fm"^ (exact 
symmetric case) and Cs = 2.0 fm"^ (non exact symmetric case). Moreover, the strength 
of the centripetal barrier term is set up to some arbitrarily chosen values: g = 0.5, 2 and 
6 corresponding to m = 0.3660254, 1 and 2, respectively. For the values g = 2 and 6, 
the singularity acts as impenetrable barrier, thus deviding the space into two independent 
regions, the negative half-line and the positive half-line. 

Dirac equation which in the limit of a non-relativistic and spinless particle transforms 
into Schrodinger equation for the isotonic potential (1) is constructed as follows. In the 
exact spin symmetry, we set Cs = and apply appropriate transformations given by 
{Mc^ + -E„,-i) /h'^c^ ~ 2M/ff and En-i — Mc^ ~ En, we finally obtain the Schrodinger 
solution in (14). 

Let us now turn to the calculations of the corresponding wave functions for this system. 
We obtain the first part 0(s) of the wave function (7) and the weight function p{s) as 

f](s) = sii+Kexp (--z/sj , (52) 
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where 



C = ^^1 + 2(77, g = m{m + 1), (53) 



and 

p{s) = s'' exp (—i^s) . (54) 

Hence, the second part yn{s) of the wave function (7) can be obtained from the weight 
function as 

y^s) ~ 4« (us) . (55) 

Finally, we find the normalized wave function satisfying Eq. (37) as 



F„,_i(x) = J ^,:^ ;':^^ xi+^exp (-^^^') 4^^ (^^') • (56) 



r(n + c + 

In addition, the corresponding lower-spinor component wave function Gn,~i{x) is found from 
the solution of Eq. (39) as 



' ^' (Mc2 + K,-i-C,)Yr(n + C + l) ^V 2 

-1 + 2C \wc). 2N , dLl^\ux') ] 

ux ] L)^' (ux) -^ ; . (57 

2x ) "^ ^ ' dx ^ ' 

Let us remark that the obtained results of the s-wave (/t = —1) of the spin-symmetric 
Dirac equation with isotonic oscillator agree with the results of the three dimensional Dirac 
equation with harmonic oscillator potential combined with the centrifugal barrier term 
K (k + 1) /r^ if we make the parameter change m -^ k since the positive half-line x G (0, oo) 
in the first is also equivalent to r G (0, oo) in the second. This is apparent because the cen- 
tripetal barrier potential m(m+l)/x^ in the isotonic oscillator is equivalent to the centrifugal 
term k{k + 1) jr^ in Eq. (37). 

On the other hand, the Klein-Gordon solution for the isotonic potential (in relativistic 
h = c = \ units) can be obtained from the exact spin-symmetric case, V = 5*, C^ = 0. 
Hence, the energy equation can be obtained from Eq. (51) as 

{El_^ - M') (E„,„i - M) = 2Mu? ( 2n + 1 + ^^1 + 2^ (Mc^ + E„,_i)j , (58) 

and the wave function from Eq. (56) as 



\ 5(1+^0) 



^"-'(^) - V r;l':'U'i) --"°-p (-r^O ^' (-^) • 
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iMa;2 (^ + ^^ _^)^ ^^ = \\l^^ '^3 {M + ^„,-i). (59) 



2. Pseudospin symmetry limit 

The exact pseudospin symmetry occurs when S ~ —V or S = Cps = constant [12,16] and 
the quahty of the pseudospin approximation in real nuclei is connected with the competition 
between the pseudo-centrifugal barrier and the pseudospin-orbital potential [32]. Therefore, 
we take the difference potential in Eq. (40) as the isotonic potential model, i.e.. 



2 2x2 



A = Ujsotix) = -MwV + -4, (60) 



In the pseudospin symmetry, the eigenstates with with j = / ± 1 are degenerate for / 7^ 0. 
For the s-wave case ( /t = 1) and in terms of the variable s = x'^, Eq. (40) reduces to a 
simple form 

G'Us) + j^fU') + ^ [^'5' - A^s + p] G„,i(s) = 0, (61) 

where 



/3 = ]^gi and V = ^J^Mu^^. (62) 

To avoid repetition in the solution of Eq. (61), a first inspection for the relationship between 
the present set of parameters (Ap^,/3,z/) and the previous set (A^,/?, z/) provides that the 
energy solution for pseudospin symmetry can be similarly found directly from those of the 
previous energy solutions for spin symmetry using the following parameters map 



U{s) ^ -U{s) (z/2 ^ -D\ p^-poiu-^ jco, g ^ -g), 3 = v^^, (63) 

from which trivial calculus gives us the transcendental energy equation: 

{En,i + Mc^) yjEn,i - Mc^ - Cps = haoV2M (2n+l + c) , (64) 

with 

C = I ^Jl + ^,{En,i-Mc^-Cps). (65) 

where n = 0, 1, 2, 3, ■ ■ ■ and En,i > Mc^ + Cps is the main condition for the real bound state 
solutions. 
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Now the normalized lower-spinor component wavefunctions of the isotonic oscillator are 
given by 



FnAx) 



\ 



2(2z/)^+^n! i^g^^^ fl~^2\ ^(c) (,^^2^ . (66) 

(n + C + l) V 2 / 



r(n + c 

Therefore, using Eq. (64), we compute some energy levels for several values of n. In the 
pseudospin symmetric limit. Table 2 gives some numerical results by taking the following 
parameters values: M = u = 1.0 fm'^, Cg = fm~^ (exact symmetric case) and Cs = —2.0 
fm^^, —13.0 fm~^ (non exact symmetric case). In addition, the strength of the centripetal 
barrier term is set up to the following arbitrarily chosen values: g = 0.5, 2, and 6. 

III. CONCLUSIONS AND OUTLOOK 

In this work, qualitative data were obtained on the modifications of spectrum energy 
on a nonrelativistic and relativistic particle confined by isotonic oscillator field of specific 
strength g. The spin and pseudospin symmetry in relativistic isotonic oscillator are inves- 
tigated systemically by solving the Dirac equation with scalar and vector radial potentials 
by applying the NU method. In one-dimensional isotonic oscillator, we have obtained the 
exact solutions in closed form for the energy spectrum and the wave functions, which are 
equivalent to solving the three-dimensional harmonic oscillator problem. The isotonic os- 
cillator is an isospectral to harmonic oscillator. Also, the energy steps are twice that of 
the simple harmonic oscillator. The resulting solutions of the wave functions are written 
in terms of the associated Laguerre polynomials L'^{z) (confluent hypergeometric functions 
M{a, b, z) = ^Fi {a; b; z)) and the wave function for states n = 0, 1 and 2 are found to have 
the same shape as the harmonic oscillator as shown in Figures 2, 3 and 4. The case where 
n is even appears to be the most interesting, since all generated wave functions are normal- 
izable. However, when n is odd, half of the generated wave functions must be removed as 
they are not normalizable. 

In the relativistic case, it is found that the solutions when A = ^ S = V {i.e., exact 
symmetric case, C^ = 0) or S = — )■ S* = —V {i.e., exact pseudosymmetric case, Cpg = 0) 
are identical to the Klein-Gordon solutions. Besides, they can be readily reduced to the 
expected nonrelativistic limit when appropriate mapping transformations of parameters are 
made. In the numerical work, the relativistic energy spectrum for the spin and pseudospin 
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symmetries are given in Tables 1 and 2, respectively. It is noticed that the parameters g, 
M, Cs and Cps should be adjusted to provide us real solutions for the energy eigenvalues 
and eigenfunctions. 

Finally, let us also mention that the isotonic oscillator possesses a remarkable property. 
The change oi x -^ ix in the wave equation resulting in the change of eigenvalues En — )■ 
—En- Then, if ipni^) is the eigenfunction corresponding to the eigenvalue En, then the 
eigenfunction ip^i^x) will be normalizable only if —En is in the point spectrum of this 
Hamiltonian. So the isotonic oscillator wave function with the change x — )■ ix would be a 
good wave function as they are normalizable (— -E„ G {n\H\n)) or would not be good if they 
are not normalizable {—En ^ {n \H\n)). The eigenfunction ipniix) can be used to generate 
new operators in the supersymmetric quantum mechanics [6,7]. This remains as an open 
question that deserves to be studied. 
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FIG. 1: Behaviour of the isotonic oscihator potential (continuous hne) and the harmonic oscihator 
potential (dash line). 



FIG. 2: Behaviour of the ground state wave function ^„^q^^^(x) of the isotonic oscillator (con- 
tinuous line) and the corresponding wave function <Pq{x) of the harmonic oscillator (dash line). 



FIG. 4: Behaviour of the second excited wave function ipj^^2m=ii^) o^ ^^^ isotonic oscillator 
(continuous line) and the corresponding wave function (^2(2^) of the harmonic oscillator (dash line). 



FIG. 3: Behaviour of the first excited wave function ipn=i m=ii^) of the isotonic oscillator (contin- 
uous line) and the corresponding wave function (piix) of the harmonic oscillator (dash line). 
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TABLE I: The spin symmetric bound state energy eigenvalues (in fm ^), for several values of n 
with parameter values M = 1.0 fm~^ and u = 1.0 fm"^. 



n/En g = 0.5 (m « 0.366) g = 2 {m = I) g = Q {m = 2) 



Cs = 2.0 /m^i 

g = 2 [m = 1) 5 = 6 (m = 2) 



2.5509860 

1 3.7292142 

2 4.7223578 

3 5.6093599 

4 6.4244044 

5 7.1861562 

6 7.9061955 

7 8.5923225 

8 9.2501029 

9 9.8836823 

10 10.4962522 



3.1503636 

4.2915849 

5.2667833 

6.1428129 

6.9503157 

7.7065008 

8.4222280 

9.1048960 

9.7598277 

10.3910117 

11.0015335 



4.0959121 

5.1735045 

6.1147629 

6.9690531 

7.7611866 

8.5058073 

9.2124501 

9.8877527 

10.5365663 

11.1625702 

11.7686371 



3.3991120 

4.6747397 

5.7095838 

6.6208542 

7.4521361 

8.2256717 

8.9547327 

9.6480343 

10.3116853 

10.9501754 

11.5669263 



4.2634174 

5.4772542 

6.4867680 

7.3835758 

8.2052891 

8.9719327 

9.6957461 

10.3848919 

11.0451537 

11.6808166 

12.2951658 



"Exact spin symmetric limit. 
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TABLE II: The pseudospin symmetric bound state energy eigenvalues (in fm~^), for several values 
of n with parameter values M = 1.0 fm~^ and to = 1.0 fm~^. 

^ps ^ Ops = —2.0 Ops = — lo.U 

n/En g = 0.5 g = 2 g = 6 9 = 0-5 g = 2 g = 6 g = 2 g = 6 

1.7353829 1.9975105 2.6220370 0.8996794 1.3991120 2.2634174 0.8228652 1.8370383 

1 2.9274128 3.2918405 3.9528022 2.1870188 2.6747397 3.4772541 1.5785297 2.5680523 

2 3.9414440 4.3370543 5.0071893 3.2260195 3.7095838 4.4867680 2.2966386 3.2659358 

3 4.8433785 5.2545579 5.9290480 4.1395244 4.6208542 5.3835758 2.9834157 3.9357442 

4 5.6693464 6.0900511 6.7671403 4.9722337 5.4521361 6.2052891 3.6435022 4.5813401 

5 6.4394382 6.8666546 7.5454937 5.7467734 6.2256717 6.9719327 4.2804724 5.2057558 

6 7.1660777 7.5980685 8.2781774 6.4765859 6.9547326 7.6957461 4.8971501 5.8114252 

7 7.8575782 8.2932428 8.9743213 7.1704749 7.6480344 8.3848919 5.4958138 6.4003383 

8 8.5198335 8.9584266 9.6402732 7.8345997 8.3116853 9.0451537 6.0783346 6.9741474 

9 9.1572079 9.5981991 10.2806717 8.4734818 8.9501754 9.6808166 6.6462725 7.5342431 

10 9.7730448 10.2160418 10.8990360 9.0905633 9.5669262 10.2951658 7.2009446 8.0818094 

"Exact pseudospin symmetric limit. 
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